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It was suggested earlier that the gravitational redshift for¬ 
mula can be invalid when the effect of the backscattering is 
strong. It is demonstrated here numerically, for an exemplary 
electromagnetic pulse that is: i) initially located very close 
to the horizon of a Schwarzschild black hole and ii) strongly 
backscattered, that a mean frequency does not obey the stan¬ 
dard redshift formula. Redshifts appear to depend on the fre¬ 
quency and there manifests a backscatter-induced blueshift in 
the outgoing radiation. 
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I. INTRODUCTION 

Standard derivations of the gravitational redshift base 
on the approximation of geometric optics (see, for a dis¬ 
cussion, Q - Q ). This problem has been recently re¬ 
considered in the wave formulation §. It was shown, 
for a class of compact shocks that are separated in some 
sense Q from the horizon of a Schwarzschild black hole, 
that the energy flux scales accordingly to the redshift for¬ 
mula. Introducing the concept of photons and assuming 
that their number is conserved, one again arrives at the 
standard relation for the frequency. As an alternative 
approach one can apply the Fourier analysis, as outlined 
below in Sec. II. 

The crucial feature behind the above mentioned com¬ 
pactness and separateness conditions is that when they 
hold true, then the backscatter is negligible [H. It was 
remarked in Q that a compact pulse of radiation that 
is exposed to a significant backscattering must not obey 
the familiar redshift relation. This is suggested by the 
following reasoning. A spatially compact wave packet 
consists of a superposition of monochromatic waves of 
different frequencies. The effect of the backscattering 
is stronger for low frequency waves than for high fre¬ 
quency ones, that is spectral amplitudes of the former 
are stronger damped than those of the latter. Therefore 
one expects that the backscattering induces a shift in 
a mean frequency, so that it does not satisfy the gravi¬ 
tational redshift formula. The aim of this investigation 
is to extend results of and find a numerical example 
in favour of this conjecture. This examplary wave pulse 
must break - and in fact does so, in accordance with 
- the compactness condition of [H. 


In the next section we write equations and review rel¬ 
evant results of the preceding papers i, i- Sec. Ill 
describes an electromagnetic wave that does not comply 
to the standard redshift formula. Sec. IV presents a 
short summary. 

II. GRAVITATIONAL REDSHIFT: CLASSICAL 
DERIVATION 

The space-time geometry is defined by the 
Schwarzschild line element, 

ds'^ =+ — ^-^dR^+ R^dn^ , ( 1 ) 

where t is a time coordinate, R is an areal radius and 
dVl^ = dO^ -b sin^ Odcf^ (with 0 < 0 < tt and 0 < (j) < 27r) 
is the line element on the unit sphere. The Newto¬ 
nian gravitational constant G and the velocity of light 
c are put equal to 1. We define the tortoise coor¬ 
dinate r* = R -\- 2mln(^ — 1) and -qn = 1 — 

We will consider only the dipole electromagnetic term 
and, more specifically, choose the potential one-form 
A = •\/3/2sin^ ,t)d(j). The unknown function 'k 

satisfies the equation ||] 

= ( 2 ) 

Let the electromagnetic strength field tensor be . The 
stress-energy tensor of the electromagnetic field reads 
— [1/A)g'^^F-^gF'^^ and the time-like trans¬ 
lational Killing vector is denoted as One can define 
the energy flux Pr, 

pR{R,t) = ^f dS{R)nrT^>^Cn- ( 3 ) 

plr JS{R) 

Here n is the unit normal to the sphere R and dS is the 
standard area element on S. Pr is equal to 

Pfl(i?,t) = (4) 

^/m 

Let L/io be a null geodesic directed outward from the 
point i?o of the initial hypersurface and let L/f^ (/{be a 
segment of T/jp that connects Rq and {R,t). Comparing 
the energy fluxes through the spheres S{R) (where R» 
2m) and the initial S{Rq), one obtains 
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PniR) « \f^pR{Ro), (5) 

V VR 

assuming that the initial support of a dipole wave 
packet is contained in the annulus (a, b) such that (b — 
a)/{ar]^) « 1. Here a and b are the areal radii. 

There are two ways to derive the standard redshift for¬ 
mula from (^). 

i) Eclectic approach. The condition 

{b-a)/{ar]l) «l (6) 


implies the validity of the geometric optics approxima¬ 
tion. Thence one can write Pr{Ro) = fiiORgNpig and 
Pr{R) = huRNR, where ur^ and a)_R are the mean fre¬ 
quencies of the initial and final pulses (as measured by 
static observers) and and Nji are the respective pho¬ 
ton fluxes. If the photon fluxes are conserved, then from 
one infers 
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( 7 ) 


ii) Classical approach. One can Fourier-analyze the 
quantity representing the strength field tensor, h = 
{—dt + this is preserved along T/ig, if condition 

is satisfied. Let the spectral strength field density 
be g{iOr) = f -I- Let the support of 

g(cu) be H. If one assumes the normalization condition 
dujTlg(‘^T)l = 1, then the average frequency can be 
defined as Cr = f^dujT‘^rlg(‘^T)l- When the analysis is 
done with respect the asymptotic time t, then asymp¬ 
totic mean frequency Cj is found. The Fourier analysis 
with respect the proper time dr = y/rf^dt of a static ob¬ 
server located at R gives a corresponding mean frequency 
ih /■ Thus the two frequencies satisfy ©■ 


III. COUNTEREXAMPLE TO THE STANDARD 
REDSHIFT FORMULA 

It was pointed out in ||^ that (|^) is valid for very com¬ 
pact initial pulses - those satisfying the assumption §)■ 
This condition demands not only that a pulse is compact, 
but also that its relative width is small in comparison to 
the relative distance from the black hole horizon. It im¬ 
plies that the approximation of the geometric optics is 
valid and that the backscatter is absent. We conjecture, 
basing on a numerical evidence, that this can be relaxed 
to {r*{b) — r*{a))/{2m) « I. 

In the case of a radiative pulse that is exposed to a sig¬ 
nificant backscattering the relation would not hold. 
Let us again review arguments in favour of this conclu¬ 
sion. A wave pulse is superposed from monochromatic 
waves of different frequencies. When resolved spectrally 
by an observer, the pulse can be seen as a collection of 
peaks, each characterized by some mean frequency. The 


effect of the backscattering is more pronounced for low 
frequency waves than for high frequency ones - the spec¬ 
tral amplitudes of the former are stronger damped than 
of the latter. Therefore one would expect that mean fre¬ 
quencies of the initial and final pulses will not 

conform to i)- The observed mean frequency Cr of an 
outgoing pulse of radiation can be blue-shifted in com¬ 
parison to the value In these circumstances 

an attempt to fit observed data of a resolved multi-peak 
spectrum to the simple scaling law of (|^) can lead to 
redshifts depending on the frequency. 

Our aim is to find numerical solutions corresponding 
to compact initial data, satisfying the compactness as¬ 
sumption b — a << a, such that: 

i) the spectral amplitude as a function of to has several 
peaks with well resolved mean frequencies; 

ii) the evolution exhibits significant backscatter (hence 
the energy flux is diminished); 

iii) some of mean frequencies at Rq and R do not com¬ 
ply to (0). 

The condition b — a « a (we remind the reader that 
a and b are areal radii) is imposed in order to guarantee 
that initially the wavelengths are much smaller than the 
radius of the wave front and the curvature radius, i.e., 
that the assumption of the geometric optics is satisfied 
at the emission region. On the other hand we do not 
assume the validity of (0); our initial data are such that 
r*{b) — r*{a) = Airm, which breaks down this inequal¬ 
ity. Let us remark that with such data one can have 
the geometric optics approximation being valid both at 
the emission and detection regions, but being evidently 
broken around R « 3m, when the wavelengths and the 
curvature radius are of the same order. 

Below we present data describing one of many anal¬ 
ysed examples. The mass is normalized to unity and 
the Schwarzschild radius Rg is equal to 2. Let I = 
{r*(a), r*(b)), the support of initial radiation, be com¬ 
prised between a = 2 -|- 1.97 x 10“^® and b = 2 + 10“^® 
or, in terms of the tortoise coordinate, between r*{b) « 
—68.46 and r*{a) « —68.46 — 47r. The purely outgoing 
initial data have the form ^ = dtfir*—t)+f{r*—t)/R 0, 
where / is a free datum in I but vanishes to the left from 
a. With I being so compact in terms of the areal radius it 
is reasonable to ignore the R - dependence in this expres¬ 
sion and to prescribe initial dt as being dependent only on 
r* — t. Thus the initial data can be defined by prescribing 

instead of /. We choose 'I'(r*) = sin^((r* — r*(a))) and 
9t'I'(r*) = — sin(2(r* — r*(a))) for r* G I = (r*(a), r*(&)) 
and 4>(r*) = = 0 outside I. 

Notice that the size of the support of the initial pulse, 
z = 47r, satisfies the inequality z > 2, which is the plau¬ 
sibility condition (according to the foregoing conjecture) 
for a significant backscatter. 

The signal is measured during the time interval 327r at 
two observation points, Tq = —68.44 and r* = 199.16; 
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in terms of the areal radius we have i?o = 2 + 1.012 x 
10“^® and R = 190.07. The grid was 120000 x 60000. 
Two snapshots of the evolving configuration taken at the 
above pair of observation points are presented in Fig. 1. 



FIG. 1. Temporal dependence of = (—dt + dr*)'^ 
(with values of h put on the ordinate) as observed at 
7?o = (2 + 10“^®)m (solid line) and R — 190.07m (broken 
line). The time (with values put on the abscissa) is measured 
from the moment of the arrival of the wave front at each of 
the observation points. The part of the radiation seen at R 
(broken line) after t = 4nm consists solely of the backscat- 
tered radiation. The signal is negligible after t > 30m - its 
amplitude is much smaller than that of the oscillatory part. 

The main pulse passes through R in less than Itt (in 
units of m) while the scattered wave appears later and 
quickly vanishes, becoming negligible after t = 30m, 
when the tail dominates. The decay of the tail asymp¬ 
totics should be like t~^, according to 0. This asymp¬ 
totic exponent was in fact obtained in many trial runs 
of our numerical code with initial pulses located outside 
a = 3m [ll| ]. 

Figure 2 in turn shows the spectral decomposition of 
h = {—dt + dr*)'^, obtained by employing the FOURCO 
package (netlib). Fast Fourier Transform (IMSL library) 
was also used, with the purpose of checking the numerical 
results coming from the application of FOURCO. We de¬ 
picture normalized frequencies in order to see clearly 
the frequency shift that is caused by the backscatter. The 
spectrum that is seen at Rq (solid line) is rescaled by 
and the spectrum seen at R (broken line) is rescaled 
by the factor In the case of negligible backscatter- 

ing both rescaled spectra would have to coincide. 

Figure 2 shows the spectral amplitudes |g(i^)| in func¬ 
tion of the frequency lo. There is a number of peaks, seen 
at both observation points. We calculated average fre¬ 
quency for each of the resolved peaks, by choosing as the 
support Q of each pulse the interval between the consecu¬ 
tive minima of |( 7 (a;)|. In the case of the lowest frequency 
peak one can observe a 25-percent blueshift, from 0.25 to 
0.31, with variances 0.1 and 0.08, respectively; the broken 
and solid lines of Fig. 2 coincide for ui > 0.5/m. In the 
remaining peaks the effect is seemingly absent. This can 
be interpreted as demonstrating that this backscatter- 
induced blueshift is frequency dependent. l/(4m) is ap¬ 


proximately the asymptotic frequency of the quasinormal 
modes. It can be seen from Fig. 2 that the effect of the 
backscattering strongly weakens modes having frequen¬ 
cies LUb (as observed by an observer static at b), such 
that LOay^ (which would be the frequency detected by 
an asymptotic static observer) is smaller than l/(8m) - a 
half of the asymptotic frequency of quasinormal modes. 



FIG. 2. Frequency amplitude |(;(ca)| (depictured on the or¬ 
dinate) as observed at Rq — {2 + 10“^®)m (solid line) and at 
R — 190.07m (broken line). The normalized frequencies 
(see the text above) are put on the x-axis (scaled in units of 
1/m). 

Modes with ujby/rjb > l/(2m) are unchanged, while 
those in the interval (l/(8m) < uJbyR% < l/(2m) un¬ 
dergo a damping of the spectral amplitude, which is the 
stronger the smaller the frequency. The physical reason 
for this is following. Our initial wave pulse has to pene¬ 
trate through the potential barrier having maximum at 
r = 3m; the modes with frequencies close to those of 
quasinormal modes are much stronger scattered and their 
transmission coefficient is small. 

Let us point out that the backscatter is negligible - due 
to the weakness of the effective potential - as long as a 
wave pulse remains in a region characterized by the tor¬ 
toise variable r* << —1. Therefore the standard redshift 
formula applies in the zone r* << 1. Define a shifted ra¬ 
diation pulse as follows: 'I'(r* + D) = d'(r*) for D < 0). 
Notice that this translated pulse will possess a differ¬ 
ent frequency profile. Let bn he defined by the equa¬ 
tion r*(6) + D = bo + 21n(6£)/(2TO) — I) and let {wt) 
be a frequency measured by a static observer located at 
r*{b). A static observer located at r*{b) + D will mea¬ 
sure a blueshifted frequency {wb)yTJbhjbD- the other 
hand a static observer located at infinity would de¬ 
tect the same set of frequencies irrespective of D. 
Moreover, this observer would notice that the deviation 
from the redshift formula (if that can be seen in the de¬ 
tected frequency range) does not depend on the position 
of the initial pulse. 

If however a physical source of a wave is fixed and 
shifted from r*{b) << —1 by some D < 0 then the 
two local sets of frequencies (u!bi,0Jbiii) measured by two 
(shifted by a distance D) local static observers do co¬ 
incide, Ubi = ujboi- (We assume that tidal effects can 
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be ignored.) The static observer at infinity would then 
see: i) two different sets of redshifted frequen¬ 
cies, (^ 62 , 00 ) and (cu 5 £, 2 ,oo); with and 

^bDi,oo = y^VbD^bDi frequencies which satisfy con¬ 

ditions ijJu,oo » 1/(4to) and Wboi.oo >> l/(4m); ii) 
those frequencies, for which uu^oo < 1/(4to) and/or 
Wfeoi.oo < 1/(4 to), would disobey the standard redshift 
formula. It is clear that for any fixed wave source with 
frequencies uji one can find a location b close to the black 
hole horizon such that some of (the would-be asymp¬ 
totically detected) frequencies are of the order of 
1/(4 to), so that (according to our numerical example) 
one could see a failure of the standard redshift formula. 

If that numerically discovered phenomen is generic, 
then there would exist a natural cutoff (of the order 
of 1/(4 to)) for the asymptotically observed frequency of 
those wave pulses that are close to an event horizon. A 
wave of a mean frequency ojb satisfying the condition 
^b^/Va « l/(4m) would not be observed by a distant 
observer. In the light of the above the heuristic analysis 
that was reported in the beginning of this section, should 
be probably supplemented by adding the following: if a 
wave peak possesses a contribution with frequencies much 
smaller than l/(4m) then (and only then) its mean fre¬ 
quency will be influenced by the backscatter. 

IV. FINAL REMARKS 

The backscattering can modify the relation between 
initial and asymptotic frequencies. A characteristic fre¬ 
quency (circa l/(4m) in the hitherto used units) that 
appears useful in this context is (in standard units) 
fc « 0.25 X 31484 x M^/m Hz, where Mq is the solar 
mass. Let us remark that fc is the frequency typical for 
the quasinormal modes of the electromagnetic radiation 
propagating in the Schwarzschild spacetime. 

A numerical example of Sec. Ill reveals a robust dif¬ 
ference between the standard redshift prediction (@) and 
the actual observation. It is observed that for a radiation 
pulse having an asymptotic frequency /oo ~ fc the initial 
frequency fu is noticeably smaller than f^ y/pR, a num¬ 
ber following from the redshift formula (^). The peaks 
with frequencies /oo » fc are seemingly not influenced 
by the backscattering. This probably is not always true 
for all radiation pulses having mean frequencies much 
bigger than fc, but we believe that in the latter case 
the deviation from the redshift formula must be small. 
Rephrasing this fact in terms of initial data, it is prob¬ 
ably fair to say that {r*{b) — r*{a) « 2m is the suf¬ 
ficient condition for the validity of the redshift formula 
( 0 )- That is a much less stringent condition than that 
used in ||^ . Let us stress that in the case of astrophysical 
black holes of stellar origin the deviation from the law 
can be seen by an asymptotic observer in the radio part 
of the electromagnetic spectrum, with wavelengths of the 


order of 100 km; we do not claim that this phenomenon 
is of astrophysical interest. 

Let us point out that if the backscattering is negligible 
then: i) all of the energy of an outgoing pulse would get 
to infinity; ii) its energy flux would be diminished, ac¬ 
cording to (Q). The only mechanism that can imply the 
loss of the radiation energy is through the backscatter 
of the radiation off the curvature of the spacetime. This 
is clearly shown in the wave formulation Q, but a con¬ 
sistent quantum field theoretic treatment (in the curved 
Schwarzschild spacetime) should give the same answer. 

In the present paper the consideration is focused only 
on the dipole term, but a similar analysis with the same 
conclusions can to be done in any multipole order. An 
analogous phenomenon, leading to the demodulation of 
gravitational signals, will also manifest in the spectra of 
gravitational waves; this in principle may be detected. 
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